Direct measurement of molecular stiffness and damping in confined water layers 
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We present direct and linear measurements of the normal stiffness and damping of a confined, 
few molecule thick water layer. The measurements were obtained by use of a small amplitude (0.36 
A), off- resonance Atomic Force Microscopy (AFM) technique. We measured stiffness and damping 
oscillations revealing up to 7 layers separated by 2.56 ± 0.20 A. Relaxation times could also be 
calculated and were found to indicate a significant slow-down of the dynamics of the system as the 
confining separation was reduced. We found that the dynamics of the system is determined not only 
by the interfacial pressure, but more significantly by solvation effects which depend on the exact 
separation of tip and surface. Thus 'solidification' seems to not be merely a result of pressure and 
confinement, but depends strongly on how commensurate the confining cavity is with the molecule 
size. We were able to model the results by starting from the simple assumption that the relaxation 
time depends linearly on the film stiffness. 



PACS numbers: 68.08.-p, 07.79.Lh, 62.10.+S, 61.30.Hn 

I. INTRODUCTION 

The structure of water, as the primary biological sol- 
vent, has been intensively studied. For example, liquid 
water adopts short-range order, which depends strongly 
on dissolved species or geometric constraints. This struc- 
ture emerges from the minimization of the free energy 
associated with the dynamic system of hydrogen bonds 
between neighboring water molecules. The entropy cost 
of the induced order almost certainly plays an important 
role in determining the structure of biological molecules 
that depend on hydration for their function, such as pro- 
teins and cell membranes^. 

A surface can act as a model system for studying these 
phenomena, as it is both geometrically disrupting and 
can be chemically functionalized to affect the structure of 
the water close to it. A particularly interesting problem is 
the emergence of density oscillations as a function of film 
thickness when water is confined between two surfaces 2 . 
This phenomenon is related to the radial density fluctua- 
tions in solvation shells of solutes. These density fluctu- 
ations have been originally observed by diffraction meth- 
ods in clay- water systems 3 ^. 

In 1982, the mechanical response of confined water lay- 
ers was directly determined using the Surface Force Ap- 
paratus (SFA)£. With the invention of the Atomic Force 
Microscopy (AFM), attempts were made to measure stiff- 
ness oscillations with this technique^. AFM probes have 
a much smaller contact area than SFA. This is an ad- 
vantage if local changes in the water structure are to be 
examined^, and potentially allows for probing regions of 
negative contact stiffness. The latter cannot be probed 
using SFA, because the instrument stiffness is not high 
enough to withstand the snap-in instability in negative 
stiffness regions. The disadvantage of AFM is that the 
signals are much smaller and the contact area is deter- 



mined by the tip shape and thus is essentially unknown. 
The small signal-to-noise ratio in AFM made the direct 
measurement of water structure an elusive goal. In 1995, 
Cleveland et al£ measured the oscillatory potential of the 
confined water layers indirectly by analyzing the Brow- 
nian noise spectrum of a AFM tip immersed in water. 
More recently, direct measurements of the structure were 
achieved by Jarvis et al£ by using nanotube probes and 
a large amplitude AFM technique, and by Antognozzi 
et al— who measured the local shear modulus using an 
AFM in shear force mode. 

In this paper we present results of direct and linear 
measurements of the normal junction stiffness of water 
confined between the AFM tip and an atomically smooth 
mica surface. This was achieved by using ultra-small am- 
plitudes of 0.36 A and sub-resonance operation, which 
avoids the problem of reduced quality factor in liquids. 
The snap-in instability was avoided by using a sufficiently 
stiff cantilever (here 0.65 N/m). This method is ideal 
to make quantitative, point-by-point measurements of 
the mechanical properties of confined water layers. The 
small amplitudes (much smaller than the nominal size of 
a water molecule) allow us to measure the elastic and 
viscous response of the confined water layer without dis- 
rupting the layers themselves, as would be the case in 
the large amplitude methods used previously. The chal- 
lenge of such a technique is the measurement of exceed- 
ingly small signals, since the usual methods of signal en- 
hancement (large amplitudes, low stiffness levers, reso- 
nance operation) are not used. Recently, we succeeded 
in implementing such a technique in UH\£iiii2*i£ - and 
in liquids 14 , using an improved fiber interferometric dis- 
placement sensor— to overcome the reduced signal-to- 
noise of the technique. Here we report on our direct 
measurements of the mechanical properties of confined 
water layers using this novel AFM technique. 



II. EXPERIMENTAL 



Small amplitude, off-resonance AFMii has recently 
been successfully used for measuring atomic bonding 



curves^*, mapping force gradients at atomic resolutio: 
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and measuring atomic scale energy dissipation^. Both 
the force gradient and the damping coefficient/ energy 
dissipation can be obtained by solving the equation of 
motion for a forced damped oscillator at a drive ampli- 
tude Aq -C A (where A is the nominal range of the mea- 
sured interaction) and w « wo. The equation of motion 
is given by: 



mx + 7± + (fcj, + k)x — k^Aa exp(icji) 



(1) 



where we linearized the force field, owing to the fact that 
lever amplitudes are much smaller than the range of the 
measured interactions, A. This assumption has recently 
been shown to be justified if Aq is sufficiently smalUt 
After solving the equation, we find for the interaction 
stiffness and the damping coefficient: 



k = kz,{—r- cos 4> — I) 



and 
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Here, Aq is the drive amplitude of the lever, A is the mea- 
sured amplitude as the surface is approached, k^ is the 
lever stiffness, <fi is the measured cantilever phase, and k 
is the measured interaction stiffness. In equation J2Jl, oj 
is the drive frequency and 7 is the damping coefficient. 
In performing the above calculations, wc intrinsically as- 
sume that elastic and viscous forces are additive, i. c. 
they act in parallel (Kelvin model). In modeling liquids, 
however, typically a Maxwell model is used in which the 
elastic and viscous (damping) term are considered to be 
in series. To convert from one to the other we can use 
the following set of equations^: 

k 2 
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where 77 and R are the viscous and elastic terms in the 
Maxwell model, respectively. 

The used AFM was home-built and incorporated a 
fiber interferometer mounted on a remote controlled 
nano-manipulator with 5 degrees of freedom and a 
step size of < 100 nm. The nominal sensitivity of the 
interferometer was 3 x 10~ 4 A/v / Hz. This allowed us 
to measure stiffnesses of a few 10~ 2 N/m using a 0.65 
N/m cantilever, a sub- Angstrom lever amplitude of 0.36 
A and reasonable integration times. The measurement 
frequency was 411 Hz. Measurements were performed 
in ultrapure water with a concentration of 0.01 M KC1. 
The surface was freshly cleaved mica, and the cantilever 
tip was made out of silicon. 




Displacement (A) 

FIG. 1: Amplitude and phase measured on a water layer con- 
fined between the AFM tip and a mica surface. The mica 
surface is located to the right. Several oscillations can be 
seen in the amplitude data. The overall decrease as the sur- 
face is approached is due to hydrophilic effects. The phase 
shows a more complicated behavior (discussed in text), but 
also shows clear oscillations. The reference lines correspond to 
displacements where liquid ordering is maximized and serve 
as a guide for the eye. 



III. RESULTS 

Figure^shows the amplitude and the cantilever phase 
as a function of displacement. The surface is located to 
the right of the graph, and the monotonic drop-off of 
the amplitude as the surface is approached can be at- 
tributed to repulsive interactions, which are most likely 
hydrophilic in origin. The amplitude data shows at least 
5 equally spaced local minima (and maxima). The phase 
data shows equally spaced maxima further away from the 
surface, which roughly line up with the minima of the 
amplitude data. However, as the surface is further ap- 
proached additional 'intermediate' peaks appear close to 
the amplitude maxima, and these peaks finally dominate 
as the gap is decreased to a few molecular spacings. The 
average spacing between the amplitude minima (and the 
phase maxima further out) is 2.56A ± 0.20A, consistent 
with earlier reportsi^^iS. Overall, the phase increases 
up to a global maximum as the surface is approached and 
then decreases again closer to the surface. 

The measured stiffness (equation (01) can be decom- 
posed into two components: A monotonic background, 
and an oscillatory term, which is the one we will be 
most concerned with in this paper. The monotonic back- 
ground is most likely due to double-layer (DVLO) and 
hydrophilic interactions, which can both be modelled as 
exponentials. A best fit and subsequent subtraction of 
the monotonic background yields the data shown in Fig- 
ure [5] Also shown in Figure |21 is the damping coefficient 
calculated using equation J2J. Note that the peaks in 
the stiffness data correspond to the minima in the ampli- 
tude data and thus to the higher stiffness of the ordered 
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FIG. 2: Solvation stiffness and (Kelvin) damping coefficient 
versus displacement. The solvation stiffness was obtained by 
calculating the stiffness from equation and subtracting the 
exponential background. Again, clear oscillations spaced at 
2.56 A can be seen in both the stiffness and the damping. The 
damping exhibits a 'phase-shift' with respect to the stiffness 
data at a displacement of 15 A. Closer to the surface the 
damping is out-of-phase with the stiffness, while further away 
it switches to being in-phase. 



phase of the confined water layer. Close to the surface, 
the damping curve shows peaks that are 'out of phase' 
with the stiffness data. Further away from the surface, 
however, double peaks occur, and finally the damping 
shows peaks that are in-phase with the stiffness maxima, 
similar to the phase data shown in Figure Due to 
the dissipation, the cantilever loses kinetic energy. The 
energy loss per cycle can be calculated fromi 9 . 



-fx = ttjujA — Trk^AoA sin < 



(G) 



cycle 



The maximum loss was -Ediss = 1-3 meV per cycle, which 
was observed close to the maximum in the phase (Figure 

EJ> 

What is the origin of the observed dissipation and its 
increase as the surface is approached? One interpreta- 
tion would be to attribute the damping to viscous drag, 
especially due to the squeezing of the liquid between the 
tip and the substrate^S. The squeeze Reynolds number, 
Re, is given by: 



Re 



pwz dz 



(7) 



where rjw is the viscosity of water, p\y is its density, and 
z is the tip-surface distance. In our case, z ~ 1 ran, 
dz/dt ^ luAk 100 nm/s and we obtain Re » 10~ 10 < 
1. In the following we assume that the viscous damping 
due to the cantilever beam does not change much with 
separation (since the beam is several micrometers away) , 
and thus the variation of the damping with separation is 
dominated by the damping at the tip. With Re<§;l, the 
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FIG. 3: Comparison of Kelvin damping, 7, and Maxwell 
damping, 77. If we treat the confined film as a liquid (Maxwell- 
type model) , we find that the damping increases very strongly 
as the film is squeezed to a few molecular layers, unlike the 
Kelvin damping which increases more moderately. The inset 
shows the Maxwell damping, 77, on a log- scale, with the stiff- 
ness as a reference. It can be seen that rj is essentially in-phase 
with the stiffness throughout the measurement range. 



squeeze damping term between tip and surface is given 
by2°: 



7 S = 6irrj w 



p2 



(8) 



Using reasonable values (10-100 nm) for the tip radius, 
^?tip, we find that the expected viscous damping at less 
than 1 nm separation is of the order of 10 -9 to 10~ 7 
Ns/m, which is about three to four orders of magnitude 
smaller than the measured values (see Figure 01, which 
are of order 10 -5 to 10 -4 Ns/m. We found that the vis- 
cosity increased exponentially with distance and thus is 
large only very close to the surface (< 1 nm). This ef- 
fect has been observed before and has been attributed 
to a sharp increase in the effective viscosity of confined 
liquid layersi2iSS*2i, a possible indication of the altered 
dynamical and structural properties of liquids under con- 
finement. However, Raviv et alm^^ recently observed 
viscosities close to the bulk value even at separations as 
small as 1 nm. Their results were based on measurements 
of the snap-in instability close to the surface, while in 
the present experiment the mechanical properties of the 
film were measured continuously without any mechani- 
cal instabilities. Moreover, in our experiment we consid- 
ered normal forces, while their results are based on shear 
measurements. How these measurements relate to each 
other and if there is a fundamental difference between 
the normal and the lateral dynamic behavior of water is 
an important question for future study. 

Another debate has been the nature of the structural 
and dynamical differences between confined and bulk liq- 
uids. The difference has been attributed to a type of first- 
order phase transformation from liquid to solid 2 ^, or, al- 
ternatively, to a continuous transformation, not unlike a 
glass transition 25 . If the liquid does indeed turn solid un- 



der certain confinement conditions, the proper mechan- 
ical model would be a Kelvin-type model. However, if 
the liquid stays essentially liquid albeit with greatly en- 
hanced viscosity (recent evidence for this comes from dif- 
fusion measurements^) , a Maxwell- type model should be 
used. As in any 'standard' analysis of AFM, we used a 
Kelvin-type model above. However, to elucidate the na- 
ture of the changes under confinement further, it is im- 
portant to use a model that properly applies to liquids. 
Using equations J3J an d © we transformed the measured 
stiffness and damping terms to the Maxwell model. We 
found that the stiffness remains almost unchanged be- 
tween the two models (i. e. k sa R), but as shown in 
Figure EH the viscous term changes dramatically. The 
Kelvin model damping term, 7, is out-of-phase with the 
stiffness oscillations close to the surface ( Figure E}, while 
the Maxwell model damping, 77, is much larger and es- 
sentially in-phase with the stiffness variations (Figured! • 
As mentioned above, further away from the surface, the 
Kelvin damping experiences a 'phase shift' and becomes 
in-phase with the stiffness (similar to the phase data in 
Figure The Maxwell damping, on the other hand, 
remains in-phase through-out. More about this below. 

When dealing with dissipative behavior it is useful to 
look at the characteristic time constants involved in the 
dynamic behavior of the confined liquid. In the Kelvin 
model, a characteristic time is given by t c — 7/fc, which 
is called the 'retardation time'—. This time is approx- 
imately the time needed to build up a significant strain 
in the material upon application of a constant stress. In 
standard solids, a certain amount of strain can be ob- 
tained almost instantaneously due to the elasticity of the 
material, however, in ideal liquids, instantaneous strain 
is not possible due to the velocity-dependent damping. 
Thus a lower t c might indicate a more solid-like material. 
On the other hand, in the Maxwell model, the character- 
istic time is t r = i]/ R, which is the 'relaxation time' . 
This time is related to the time needed for stresses in 
the material to relax after a strain has been imposed. In 
solids, stresses will persist for long times when a strain 
is applied (one of the characteristic of materials being 
solid), while in liquids any stresses will quickly dissipate 
away. Thus higher t r indicates a more solid-like behav- 
ior. It should be noted that t c and t r are simply related 
by: 

t c = \t r (9) 

The dependence of t c and t r on displacement is shown in 
Figure It can be seen that overall t c is decreasing and 
t r increasing as the liquid layer is increasingly confined. 
This indicates a tendency for the layer to become more 
solid-like. Even more interesting, however, is the fact 
that at separations where the stiffness oscillations are at 
their maximum, t c is lowered and t r is increased. This is 
a further indication that the water becomes more solid- 
like when it is allowed to order, i. e. when the tip-surface 
separation is commensurate with the 'natural' molecular 
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FIG. 4: Relaxation (t r ) and retardation (t c ) time plotted ver- 
sus displacement. The relaxation time is in-phase with the 
stiffness, while the retardation time is out-of-phase. The re- 
laxation time increases overall due to the increasing pressure, 
but also exhibits in-phase oscillations due to solvation effects. 

spacing of water. In the 'ordered phase' , the stiffness k 
(or R) is maximum, the retardation time, t c , is minimum, 
and the relaxation time, t r , is maximum, as expected for 
a solid. 

Another interesting observation is that the character- 
istic times t r and t c change very slowly with separation 
until a few Angstrom from the surface. Then they 
seem to change more rapidly, with the liquid becoming 
seemingly even more solid-like. Some authors^* have 
argued that this behavior shows the liquid undergoes 
some kind of first order phase transformation upon 
confinement, while others suggest a more gradual, glass- 
like transition— Recent results measured under shear 
even suggest that water in particular fails to 'solidify' 
at all2i. Based on our results we cannot conclusively 
decide between these viewpoints. However, as we will 
see below, there seems to be both a gradual stiffening 
of the layer as pressure is applied and a much more 
pronounced periodic change of the mechanical behavior 
of the layer. 



IV. DISCUSSION AND MODEL 

There are several surprising findings from the linear 
measurement of confined water presented here: 1) Ob- 
served oscillations in the phase and dissipation extend 
much further than the oscillations in the amplitude or 
stiffness, 2) the phase and the Kelvin damping oscilla- 
tions experience a 'phase shift' with respect to the stiff- 
ness data as we move away from the surface, 3) while on 
average the amplitude continuously decreases as the sur- 
face is approached, the phase seems to pass through an 
intermediate maximum, and, finally, as hinted above, 4) 
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the mechanical behavior of the layer changes both gradu- 
ally (hydrophilic background) and more abruptly (solva- 
tion shell oscillations). To explain this behavior we simu- 
lated the nanomechanical behavior of the water layer by 
starting from the assumption that the relaxation time, 
t r , is to first order linearly dependent on the stiffness of 
the water layer. This is not to be taken literally, in the 
sense of a direct physical connection between the stiff- 
ness and the relaxation time (although there well might 
be), but rather the stiffness is seen as an indicator of the 
'solidncss' of the layer, and the relaxation time (as an- 
other indicator) is taken to be essentially proportional to 
it. We found that we can get the best fit of our data if 
we assume that the relaxation time depends linearly on 
both the background stiffness, kh, due to hydrophilic in- 
teraction and on the stiffness oscillations, fc s , due to the 
solvation effects but with two different 'coupling' con- 
stants ai and cti'- 

t r = u\k h + a 2 k s + t a (10) 

Here, to corresponds to the relaxation time measured far 
away from the surface. The advantage of using separate 
constants ot\ and cti is that we can separate the effect of 
background hydrophilic interactions from the influence of 
solvation forces on the relaxation time. We found that in 
order to reproduce the experimental results as closely as 
possible it was necessary to set a\ to 2.3 x 10 -3 sm/N and 
Oi2 to 5 x 1CP 3 sm/N, i. e. the relaxation time was more 
than twice as sensitive to solvation forces than it was to 
the hydrophilic background. It should be noted that the 
hydrophilic background stiffness is directly proportional 
to the load (or surface pressure), since both are exponen- 
tials and one is the derivative of the other. Thus by the 
above approach we can separate the effects of the over- 
all pressure or load from the effect of the liquid ordering 
which only occurs at certain, molecularly commensurate 
separations. 

The retardation time, t c , can be calculated from equa- 
tion . The damping coefficients 7 and r\ are then given 
by: 

l = t c -k (11) 
r) = t r -R (12) 

In the simulations, we took k — R (as found experimen- 
tally) to simplify the calculations. The solvation force 
was modeled as follows 2 : 

N 

F s = ^r 2 {z k )- 

tip,fc=0 

{ 2n(z k + D) \ ( z k + D \ 
■k B Tp cos I I exp I I (13) 

where we summed the contributions of different areas of 
the tip by subdividing the tip into N horizontal 'slices' 
of radius r(z k ). Here, p is the particle density of water, 
a is both the period and the decay parameter of the in- 
teraction (they were experimentally found to be nearly 



identical), z k is the height of the fc'th slice of the tip, and 
D is the tip-surface separation. The hydrophilic interac- 
tion is given by: 

F h = J2 2nr 2 {z k )-p h exp(-^^-j (14) 

tip,fc=0 ^ ' 

where pn is a constant and A is the decay parameter of the 
hydrophilic background. All parameters in expressions 
(|13fl and l|14fl were determined from the experiment. The 
hydrophilic decay parameter was found to be only slightly 
smaller (A = 2.45A) than the decay parameter of the 
oscillations (a = 2. 56 A). The corresponding stiffnesses 
were found from taking the derivative of the forces with 
respect to tip-surface separation, k = — dF/dD. 

Since we cannot know the exact geometry of the tip, 
we did not expect to get a perfect agreement between 
theory and experiment. Nevertheless we obtained a semi- 
quantitative agreement that reproduces all of the surpris- 
ing features mentioned above. The geometry of the tip 
was assumed to be paraboloid, and the best agreement 
with experimental data was obtained for a nominal radius 
of lnm. Figure |5l shows the calculated total stiffness and 
Kelvin damping coefficient, 7. The damping coefficient is 
out-of-phase with the stiffness close to the surface, then 
undergoes a 'phase shift ' and becomes in-phase close to 
the surface, as seen in the experiment. From equation 
(|llf) we see that the damping is a product of the retarda- 
tion time, t c and the stiffness, k. Close to the surface the 
damping is dominated by the retardation time t c , which 
is always out-of-phase with the stiffness (Figure while 
further away from the surface the stiffness k dominates 
the variation in the damping. On the other hand, the 
Maxwell damping, r), is always in phase with the stiff- 
ness, since the relaxation time t r is in-phase with R (or 
k). 

It can be shown that a more complicated mechani- 
cal model, such as the commonly used Burger's model, 
behaves like a Maxwell model at low frequencies. This 
implies that the present discussion has more general im- 
plications than might be expected from the use of such 
simplified mechanical models. In particular, it would 
seem from equation <|llf) that the oscillatory behavior 
of the Kelvin damping, 7, could be explained by the os- 
cillatory behavior of the stiffness, even if the retardation 
time is constant or slowly varying. In this scenario, the 
observed oscillations of the retardation/relaxation time 
would be merely an 'artifact' of the calculation. How- 
ever, if the retardation time were constant or smoothly 
varying (i.e. not oscillating), the Kelvin damping would 
have to remain in-phase with the stiffness at all times. 
This is not observed in the experiment. The fact that 7 
is out-of-phase close to the surface implies that the relax- 
ation/retardation time of the liquid exhibits separation- 
dependent oscillations independent of the oscillations of 
the stiffness. This means that the dynamic behavior of 
the confined liquid is strongly affected by how commen- 
surate the tip-surface spacing is with regard to the size 
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FIG. 5: Simulated stiffness and Kelvin damping, 7. Compare 
to measured data (Figure Although we had to assume 
larger stiffness oscillations in the model, the overall agree- 
ment is good, and the 'phase shift' in the damping data is 
reproduced well. 
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FIG. 6: Simulated amplitude and phase. Compare to mea- 
sured data (Figure 0. There is good qualitative agreement 
and the complicated behavior of the phase is well reproduced 
including the 'phase shift' with respect to the stiffness (here: 
amplitude) and the global maximum. 



of the confined molecules. 

The phase was calculated by solving equations |J2J and 
© simultaneously (and assuming uj -C ojq): 

tan <f> = - U1 l (15) 
k + Kl 

The simulated phase, 0, is shown in Figure The sim- 
ulation reproduces all the 'puzzling' features of the ex- 
periment: The shift from being out-of-phase to being in- 
phase with the stiffness oscillations and the intermediate 
maximum in the phase. The shift is due to the shift in 
7 discussed above. The intermediate maximum is due to 
the fact that the stiffness changes slowly far from the sur- 
face, but then rather rapidly closer in, 'overtaking' the 
damping coefficient in the process (equation j!5|>). The 
observation that oscillations in the phase or damping are 
observable further away from the surface than the oscil- 
lations of the stiffness can also be explained: As we can 
see in the simulated phase, a phase of more than 1° is 
observed as far away as 13A from the closest approach 
(about 5 water layers). Such a phase angle can be easily 
measured with a lock-in amplifier. On the other hand, 
at the same separation, the stiffness is only 0.04 N/m 
requiring a measurement of a change in lever amplitude 
of the order of 0.02A, which is more difficult to measure 
and can be lost in the noise. 

In conclusion, we can see that our simple approach of 
directly relating the relaxation time to the stiffness of 
the layer has allowed us to reproduce all the important 
features of the experiment. The relaxation/retardation 
times can therefore be taken as fundamental physical 
parameters (together with the stiffness) that character- 
ize the mechanical properties of the system quite well. 
The weaker dependence of the relaxation time on the 
hydrophilic interaction and the more pronounced depen- 
dence on solvation forces, suggests a compromise in the 
continuing debate over the nature of the solid- liquid tran- 
sition. It seems that there is a gradual increase of the 
relaxation time with surface pressure and a more sub- 
stantial change related to the molecular ordering of the 
liquid close to the surface. 
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